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Abstract 

Perfect simulation of an one-dimensional loss network on R with length distribution 7r and 

> : 

cable capacity C is performed using the clan of ancestors method. Domination of the clan of 
ancestors by a branching process with longer memory improves the sufficient conditions for the 
^ ■ perfect scheme to be applicable. 
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1 Introduction 

Kelly (1991) introduced a continuous unbounded loss network described as follows. Imagine that 
users are arranged along an infinitely long cable and that a call between two points on the cable 
si, S2 € R involves just that section of the cable between s± and S2- Past any point along its 
length the cable has the capacity to carry simultaneously up to C calls: a call attempt between 
s\, S2 G R, si < S2, is lost if past any point of the interval [s±, S2] the cable is already carrying C 
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calls. Suppose that calls are attempted at points in R following a homogeneous Poisson process 
with rate A. Assume that the section of the cable demanded by a call has distribution tt with 
finite mean p\ and the duration of a call has exponential distribution with mean one. Assume 
that the location of a call, the cable section needed and its duration are independent. Let m{s, t) 
be the number of calls in progress past point s on the cable at time t. Kelly (1991) conjectured 
that ((m(s,t),s € R),i > 0) has a unique invariant measure, given by a stationary M/G/oo 
queue (Markov arrivals, general service time and infinite servers) conditioned to have at most 
C clients at all times. Ferrari and Garcia (1998) used a continuous (non-oriented) percolation 
argument to prove the above conjecture whenever tt has finite third moment and the arrival 
rate A is sufficiently small. Fernandez, Ferrari and Garcia (2002) using an oriented percolation 
argument improved this bound to 

x ffg _ 1 ci i \ 

where p\ and p2 are the first and second moment of distribution ir respectively. This argument 
is based on a graphical representation of the birth and death process and it is the basis for the 
perfect simulation scheme 11 Backward-Forward Algorithm" , described in Fernandez, Ferrari and 
Garcia (2002). This algorithm involves the "thinning" of a marked Poisson process — the free 
process — which dominates the birth-and-death process, and it involves a time-backward and a 
time- forward sweep. The initial stage of the construction is done toward the past, starting with 
a finite window and retrospectively looking to ancestors, namely to those births in the past that 
could have (had) an influence on the current birth. The construction of the clan of ancestors 
constitutes the time-backward sweep of the algorithm. Once this clan is completely constructed, 
the algorithm proceeds in a time-forward fashion "cleaning up" successive generations according 
to appropriate penalization schemes. The relation "being ancestor of" induces a backward in 
time contact/ oriented percolation process. The algorithm is applicable as long as this oriented 
percolation process is sub-critical. Garcia and Marie (2003) using the Perron-Frobenius theory 
for sub-criticality of branching process obtained a new bound given by 

K = T-^ 1 7- (1-2) 

However, studying the characteristics of the clan of ancestors through simulation in Section 5 
it is clear that the domination by the branching process is not sharp. That is, the number of 



ancestors is much smaller than the total number of the population in the branching process 
and the clan of ancestors can be finite even though the branching is supercritical. By studying 
the perfect simulation algorithm and constructing a 2-generation branching process, that is a 
Markov process with order 2, it is possible to improve bound (1.2) to 



and it stands that A** > (4/3) A* . 

2 Loss networks 

Let Q be a family of intervals of the line 7 (7 = (x,x + u),x,u G R), which will be named calls, 
and consider a state space S = {£ G N e : £(7) / only for a countable set of 7 G Q}. 

Loss networks are a particular class of spatial birth-and-death processes. The evolution of 
these processes in time are given either by the birth of a new call to be added to the actual 
configuration or by the death of an existing call that will be eliminated from the actual config- 
uration. Moreover, they have the Markovian property in time that, the probability of a change 
depends only on the actual configuration of the system. A loss network rjt is defined by a marked 
Poisson process where births are controlled by a birth rate, a non-negative measurable function 
b(7, 77) such that 



for each B, bounded Borel set, and for all rj G S. The births are regulated by the exclusion 
principle, depending on the capacity (C) of the network. The marks include a life-time expo- 
nentially distributed with mean one and the length of the call which is distributed according 
to a distribution it with finite first and second moments p\ and P2 respectively. The death rate 
also a non-negative measurable function d : M. d x S — > [0, 00). In this work, we are going to 
assume that the death rate is always equal to one. The generator of the process is given by 




4 



(1.3) 



3/oi + VpI + 8 P2 





where r/ G {0, The death rate 1 is included in the second expression. 
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For a process at with rate densities which are independent of the actual configuration there 
exists to : Q — > [0, oo) such that 

b( 7 ,a)=^( 7 ) (2.2) 

we call this process a free process. Such a process is just a space-time marked Poisson process. 
It exists and is ergodic whichever the choice of w. In the particular case where a; (7) = A the 
invariant measure is the A-homogeneous Poisson process. For the one-dimensional loss networks 
(described in Section 1) the birth rate is uniformly bounded and it can be decomposed as 

b((x, x + u), rj) = A vr(u) M((x, x + u),rj) (2.3) 

where, < M(j,r)) < 1. The first factor represents a basic birth-rate density due to an 
"internal" Poissonian clock and the last factor acts as an unnormalized probability for the 
individual to be actually born once the internal clock has rang. The birth is hindered or 
reinforced according to the configuration rj. 
In this case, for capacity C = 1, 

M( 7 ,r?)= [J (i-AT^)) (2-4) 



1(7,0) = { 



1 7 n#/0 

otherwise 

where 7, 6 are of the form (x, x + u). For C > 1, the expression is less simple 



(2.5) 



1 otherwise 



M((x,x + u),rj) = < 



there exists y £ (x, x + u) and 
9i,... ,0c such that ij(0i) = 1 
and y £ Oi for all i = 1, ... , C. 



3 Graphical construction for the loss networks 

Let N = { (£1, Ti), (£2, 12),. . . } be a homogeneous Poisson Process with rate A in R x [0, 00), 
Si, S2, ■ ■ ■ be i.i.d. random variables exponentially distributed with mean one and U\, U2, ■ ■ ■ be 
i.i.d. random variables with common distribution it. Assume the family of variables {Si, S2, ■ ■ •}, 
{Ui, U2, ■ ■ •} and the Poisson process are all independent. Consider the random rectangles 
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Ri = {{x,y);& < x < & + u^Ti < y < r» + Si}. 

Then {Ri, i > 1} = {(&, 1}) + A, i > 1} is a Boolean model in R 2 where A = [0, x [0, Si] 
and represents the free process of calls. Boolean models have the property that the number of 
sets C G C that cover a fixed point x G M. d is a Poisson random variable with mean AE(vol(5)). 
For more details about coverage processes see Hall (1988). 

Now, for each rectangle Ri we associate an independent mark Zi ~ U(0, 1), and each marked 
rectangle we identify with the marked point Tj, Si,Ui, Zi). We recognize in the marked point 
process R = Tj, Si, Ui, Zi), i = 1, 2, . . .} a graphical representation of the birth and death 
process with constant birth rate A, and constant death rate, equal to 1. We call this free process 
a and Z; L will serve as a flag of allowed births. Calling R = (£, r, s, u, z), we use the notation 

Basis (R) = (£,£ + u), Birth (i?) = r, Life (i?) = [r, r + s], Flag (J?) = z. (3.1) 

We also define, for two rectangles R and i?', 

# * R, if # n i? + 
i?' ~ i?, otherwise. 

We need a series of definitions: 

• For an arbitrary point (x, t) G R 2 define the collection of all rectangles in R that contain 
this point 

A (x,t) = |£ G R | x G Bagis ^ 4 G Life ^| ^ 

• For each rectangle i? define its ancestor set 

Af = {R 1 G R| Birth (#) < Birth (R), R' *> R} (3.3) 

• Define recursively the generations (n > 1) of the above sets that is, the nth generation of 
ancestors: 

A**'*) = {R"\R" G Af' for some R' G A<fj$} (3.4) 
Af = {R"\R" G Af' for some R' G Af_J (3.5) 

We say that there is backward oriented percolation if there exists one point (x, t) such that 

(x t} 

An 7^ for all n, that is, if there exists one point with an infinite number of ancestors. 
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Call clan of ancestors of (x, t) the union of all its ancestors: 

A (x,t) = J A (x,t) ^ 

and R[0,t] = {R G R| Birth (R) G [0,t]}. 

The existence of the process in infinite volume for any time interval is guaranteed as long as 
the process do not explode, that is, no rectangle has an infinite number of ancestors in a finite 
time. The following theorem is proved in Fernandez, Ferrari and Garcia (2001). 

Theorem 3.7 If A^'*) n R[0, t] is finite with probability one, for any i£R and t > 0, then for 
all A C R the loss network process defined in A is well-defined and has at least one invariant 
measure /i A . 

For the existence of the process in infinite time, it is needed that the clan of ancestors of all 
rectangles are finite with probability one, that is, there is no backward oriented percolation. In 
order to construct the invariant measure for stationary Markov processes it is usual to construct 
the process beginning at — oo with an arbitrary configuration and look at the process at time 
0. If the configuration at time does not depend on the initial configuration then we have a 
sample of invariant measure. The graphical construction described above allow us to construct 
the process rjt by a thinning of the free process at for all t G R. Moreover, the same argument 
shows that the distribution of r/o does not depend on the initial configuration. The next theorem 
summarizes the results about the process, see Fernandez et al. (2001, 2002). 

Theorem 3.8 If with probability one there is no backward oriented percolation in H, then the 
loss network process can be constructed in (— oo, oo) in such a way that the marginal distribution 
of rjt is invariant. Moreover, this distribution is unique and the velocity of convergence is 
exponential. 

One way of determining the lack of percolation is the domination through a branching 
process. Establishing sub-criticality conditions for the branching process we obtain sufficient 
conditions for lack of percolation. Looking backward, the ancestors will be the branches. The 
time of the death will be the birth time for the branching process. The clan of ancestors in 
itself is not a branching process because the lack of independence. 
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Let R be a rectangle with basis 7 = (x,x + u) with length u, born at time 0. Define b^(v) 
as the number of rectangles in the nth generation of ancestors of R having basis with length v: 

6» = \{R' € A*| |Basis(fl')|=t;}|. (3.9) 

The process b n is not a Galton- Watson process but it can be dominated by one (call it b n ) 
as described by Fernandez et al. (2001), where each call length represents a type and it has 
as offspring distribution the same one as b\. The number of types can be finite, countable or 
uncountable depending upon the distribution tt. 

For the one-dimensional loss network the offspring distribution of b n is Poisson distributed 
with mean 

m(u, v) = A tt(v) (u + v) (3.10) 

where m(u, v) is the mean number of children type v for parents type u. In this case, Garcia and 
Marie (2003) used the Perron-Frobenius theory for sub-criticality of branching process obtained 
a new bound given by 

A(V^ + Pi)<l. (3-11) 

In Section 6 we are going to improve this bound dominating the clan of ancestors by a 
branching process with order 2, that is, the reproducing mechanism is governed not only by the 
parents but also by the grandparents. 

4 Backward- Forward Algorithm (BFA) applied to loss 
networks 

The BFA was introduced by Fernandez, Ferrari e Garcia (2002) to perfect simulate from spatial 
point processes which are absolutely continuous with respect to a Poisson point process and 
that are invariant measures of spatial birth and death processes. 

The algorithm does involve the "thinning" of a marked Poisson process — the free process— 
which dominates the birth-and-death process, and it involves a time-backward and a time- 
forward sweep. But these procedures are performed in a form quite different from previous 
algorithms. The initial stage of our construction is done toward the past, starting with a finite 
window and retrospectively looking to ancestors, namely to those births in the past that could 
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have (had) an influence on the current birth. The construction of the clan of ancestors consti- 
tutes the time-backward sweep of the algorithm. Once this clan is completely constructed, the 
algorithm proceeds in a time-forward fashion "cleaning up" successive generations according to 
appropriate penalization schemes. 

The relation "being ancestor of" induces a backward in time contact/oriented percolation 
process. The algorithm is applicable as long as this oriented percolation process is sub-critical. 

To simplify the implementation of BFA to the loss network process we are going to assume 
that 7r have compact support. This assumption is not necessary and can be removed with a 
little modification on the generation of the free process. Define H = inf{y | 7r((0, y)) = 1}. 

4.1 Construction of the clan of ancestors of a finite window 

A = [a, b] C R 

We are interested in sampling a finite window A = [a, b] of the equilibrium measure in infinite- 
volume. 

CI. Generate the free process «o = {£1 , • • • , £m}; a homogeneous Poisson process with rate A 
in the interval [a — H,b]. 

S L = °; s R = b - 

C2. Generate Uf, . . . , i.i.d. random variables with common distribution n and let r\ = 0. 
For each i from 1 to m 

if (4°, 4° + V?) n [a, b] + then r, = r, U £° + U?) (4.1) 

We are simply generating rectangles with basis intersecting [a, b]. We have no = \rj\ < m 
basis. 

C3. Generate S®, . . . , S® i.i.d. exponential random variables with mean one and construct 
the rectangles 

Ro={fe°, $ + U?) x 0]; i = 1, . . . , n}. (4.2) 

Consider now the following subset of E x (— oo, 0] 

n 

A = {J($-H,$ + U?)x[-S?,0} (4.3) 
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C4. k = 1; A = A ; 



C5. s k L = mm(s k L \mm i < nk _ 1 (^ 1 - H)) 

4 = m a x(4" 1 ,max l < rife _ 1 (4 fe - 1 + U^ 1 )) 

C6. Generate a A-homogeneous Poisson process , T k ), . . . , , T^ fc )} on A U [s|, s^" 1 ) U 

C7. Generate C/f , . . . , i.i.d. random variables with distribution ir and S k , . . . , S k k i.i.d. 
exponential random variables with mean one and construct the rectangles 

R*={(tf,tf + Etf) x [t? - S?,r k ];i = 1, . . . ,n fc }. (4-4) 

Consider 

A fe = - H > e* + ^) x - ^ °]- ( 4 - 5 ) 

1=1 

C8. - if rife = then construct the clan of ancestors of r\ 

k-l 

A" := (J Hi (4.6) 

i=0 

and STOP. 

- otherwise, do 
A = A fc \A fc _u 
k=k+l; 
return to C5; 

Remark. At Step C6., it is necessary to consider only rectangles satisfying 

rf-S k > . min (rf" 1 - Sf- 1 ). 

j=l,...,n fc _i 

This restriction does not affect the generation of the clan of ancestors but reduces drastically 
the computational cost. 

We finish performing the BACKWARD step of the algorithm: the construction of the 
clan of ancestors. The FORWARD step corresponds to move to the beginning of the clan of 
ancestors and decide which rectangles are going to be kept and which ones are going to be 
erased. Once these clans are perfectly simulated, it is only necessary to apply the deterministic 
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"cleaning procedure" , based on the capacity C of the network , to obtain a perfect sample of 
the interacting process. In this case, if a point (x, t) belongs to more than C rectangles, keep 
the C rectangles born first and erase the others. 



4.2 The cleaning algorithm 

Call T the set of rectangles to be tested and K the set of kept rectangles. 

LI. K=0; T= A T >; 

L2. If T = go to L4. 

otherwise, order T by birth time. Let R\ be the first rectangle following such ordering. 
K = KURi; T = T\Ri 

L3. Depending upon C 

1. If C = 1; For all R G T such that R oo R 1 , T = T\R. 
return to L2. 

2. Se C > 1; 

for i=l to |T| - C 

Ri G T, if Ri oo R 1 call Area=R, t n i?i, C(Area)=2 and K = K U R; T = T\R; 
/or j=7 to |T| 

if ^ n Area / take C(Area) = C(Area) + 1, if C(Area) > C then T = T\Rj 
return to L2. 

L4. Take = K and STOP. 



Obtaining K^, we define 

( 7 ) = < 



1 7/(7) = 1 and 3R G K T ' such that Basis (R) = 7 
otherwise 



(4.7) 



Theorem 3.18 of Fernandez et al. (2002) guarantees that rf is a perfect sample from the invariant 
measure of the loss network described above. 
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4.3 Simulation results 

In this section we present some of the simulation results for several values of A, C ( network 
capacity) and window A. The distribution it is taken to be U(0, 1). In this case, Garcia and 
Marie (2003) obtained that A < 0.9282 is a sufficient condition for the simulation. The programs 
were written in MATLAB 5.0. For easiness of reading the results are presented in two steps: 
the clan of ancestors and the cleaning result. The basis of the rectangles kept at time t = 
constitutes the perfect sample. 

12 | 1 1 I =; 1 1 1 1 1 



10 - 



8 - 



6 - 



4 - 



2 - 



I U 1 1 U 1 1 U 1 — I 

123456789 

Figure 4.1: Clan of ancestors for U(0, 1), A = 0.7, A = [0, 10]. 
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Figure 4.2: Cleaning procedure for the clan presented in Figure 4.3. a) C — 1 b) C — 2. 
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5 Studying the characteristics of the clan of ances- 
tors through simulation results 

We perform a 1,000 simulations for several values of A < A*. Conditioned on the event "the 
point (x, 0) is present at the free process" (which has probability 1 — exp{— Ap}), we observed 
the values of N(A^ X ' ^) - total number of rectangles present in the clan. 

The expectations of this variable were estimated through the sample mean and compared 
then to the expected values for the associated branching process used to find the sub-criticality 
condition. 

The simulations were performed in two cases, when ir is the U (0, 1) distribution and when 
7r is concentrated in one point (fixed call length). From Garcia and Marie (2003), the critical 
value for A to assure sub-criticality is 

• When 7T = [7(0, 1) 

A* = — = « 0.9282 (5.1) 

i i /T 

• When ir(d) = 1 



2 T V3 



K = y d (5.2) 

Graphicl: N (uniform) 

Figure 5 shows that the branching process dominates the clan of ancestors (we constructed 
then this way). However, due to the fact that in the branching process we can have subsequent 
generations of rectangles to be born in the same area as the predecessor generations, the number 
of rectangles of the branching process grows much faster than the number of rectangles of the 
clan of ancestors, as A increases. 

5.1 Estimation of the critical value using simulations 

The purpose of this section is to study the behavior of the clan of ancestors as A increases above 
A*. From Figure 5 we can see that the critical value obtained through the domination by a 
branching process underestimates the true value of the finiteness of the clan of ancestors. The 
idea behind these results is to generate samples for increasing values of A and to study the total 
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Branching process 
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Clan of ancestors 




Branching process 



» 8 8 8 S ooooooooc 



0.0 



0.2 



0.4 0.6 
Lambda 



0.8 



1.0 



Figure 5.3: Expected total number of rectangles in the branching process and the clan of ancestors, 
a) 17(0,1) b) d = 0.5 
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number of rectangles. Our conjecture is that, close to the true critical value A c the total number 
of rectangles should grow exponentially fast. Thus finding an assintote for E(iV) would give us 
an estimate of A c . This is true for the branching process, comparing the value of E(iV) as A 
approaches A* in Figure 5 we can see visually a vertical assintote at A*. 

From now on, for all distributions, we sampled 1,000 observations of the clan of ancestors 
and computed N (the sample mean) for each value of A. Figure 5.1 present the results for fixed 
length calls, d = 0.5. 



600 - 

z arc 

400 - 

- 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Lambda 

Figure 5.4: Expected number of rectangles in the clan of ancestors (E(iV)) for d = 0.5 

At first sight we see that there is an assintote close to 2.8. To be more precise, we tried to 
find a root for the equation l/log(7V) = . We used a degree 19 polynomial to approximate 
1/ log (iV) and found a root in A = 2.8231. 

Several simulations were performed for several values of d just to get a more precise estimate 
for A c since due to the invariance of the Poisson process for fixed call length there is a linear 
relationship among the critical values for all d, see Table 5.1. 

Comparing the values of A* = ^ and A c we can see, as expected, a linear tendency and we 
can adjust a regression model with no intercept using least squares to get 

A c = 2.8246 (5.3) 



□ooooc 
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d 


0.5 


1.0 


1.4 


2.0 


2.5 


3.0 


3.5 


4.0 


4.5 


5.0 


A c 


2.8231 


1.4193 


1.0254 


0.7312 


0.5682 


0.4537 


0.3931 


0.3530 


0.3103 


0.2833 



Table 5.1: Critical value of A c obtained through simulation for several call lengths 



The question now is to perform the same comparison using different random distributions for 
7r. We simulated clan of ancestors for several Beta distributions and compared and A c . Table 
5.1 presents these results along with the ratio A c /A*. We can see that A c « 2.82A* . Adjusting 
a least square model without any intercept: 

A c = 2.8243 A: . (5.4) 



Distribution 


A c 


Ac/a: 


U(0,1) 


2.6135 


2.8157 


Beta(2,l) 


2.0888 


2.8695 


Beta(2,2) 


2.6746 


2.8022 


Beta(3,l) 


1.8597 


2.8353 


Beta(3,2) 


2.3079 


2.8444 


Beta(l,2) 


3.7981 


2.8166 



Table 5.2: Critical value of A c obtained through simulation for distributions U(0, 1) and Beta(a,/3) 



6 Dominating the clan of ancestors by a 2-generation 
branching process. Critical value. 

To construct the Galton- Watson process b n mentioned by the end of Section 3, Fernandez 
et al. (2001) used a multitype branching process B n , in the set of cylinders, which dominates 
A n . To do this they looked "backward in time" and let "ancestors" play the role of "branches". 
In particular, births in the original marked Poisson process correspond to disappearance of 
branches. Like them, we reserve the words "birth" and "death" for the original forward-time 
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Poisson process. This construction can be done by enlarging the probability space and defining, 
for any given set {Ri, ■ ■ ■ , Rk}, independent random sets B^ 1 with the same marginal distribu- 
tion as A 1 1 . The important point here is that 

k k 

[jAf*c[jBf l . (6.1) 

i=i i=i 

The procedure defined by Bi naturally induces a multitype branching process in the space 
of rectangles. The n-th generation of the branching process is defined by 

B* = {Bf : R' e B*_i} (6.2) 

where for all R', Bf' has the same distribution as Af' and are independent random sets de- 
pending only on R'. Then, for all n > 1 it follows 

A« C B* (6.3) 

It is introduced then a multitype branching process in the set of calls Q. For an initial 
rectangle R with the basis of size u, b^ is defined as the number of rectangles in the n-th 
generation of B„ that have basis with length v: 

b» = \{R' € B*| |Basis(i2')|=«}|. (6.4) 

The following relation is then established 

|A*| < |B*| = (6.5) 

V 

We will have to further enlarge the set of types of B„ by creating a branching process in the 
set of cylinders B* with a new type called Color G {b, g} which is inherited by a cylinder not 
only from the types of the first generation (parent) but also from the 2nd generation (grand- 
parent). The types b,g are chosen to associate to black and green. 

The idea is to identify some rectangles in the branching process B„ that could never be present 
in the A n . These rectangles we have in mind as "black". Then, if the total number of "green" 
(not-black) rectangles is finite implies that the clan of ancestors is also finite. We expect then 
that counting only "green" rectangles one can obtain better critical value for A. That is what 
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we are going to prove in the following. 



For a rectangle V = (£y, ry, sy, «y) in the nth generation of the branching process B„, we 
define an extra type to be its color as follows. For n = 1 or 2, Color (V) = g. For n > 2, 
let R G B*_ 1 and i?' G B*_ 2 such that V G Bf and i? G Bf' . Denote i? = (£,t,s,u,c) and 
i?' = (£', r ', s', u', d), where c = Color (R) and c' = Color (R f ). Therefore, 

• if c = b then Color (V) = b, 

• if c = g then Color (V) = 6 if and only if (£y,ry + sy) G D(uy, R, R') where 
D(u v ,R,R') = L(£,u,£,u',u v ) x [0,r'] and 

L(£, u, u', u v ) = [max(( - u v ,£' - u v ), min(£ + u, ^' + u')]. (6.6) 

This way, for every rectangle is defined its bi-dimensional type Type(-) G Q x {b, g}. 

Consider d^" ,c \(u, c), (v, c")) to be a 2-generation multitype branching process defined by 

dt c '\(u,c),(v,d')) (6.7) 
= \{V : V G B% and i? G Type (C) = («;, c'), Type (i?) = (n, c), Type (V) = (v, d')}\, 

that is, it is the number of cylinders of type (v, c") in the nth generation of a cylinder of 
type (w,d) which have an ancestor in the n — 1 generation of type (u,c). We can think of 
this branching process as a first order branching process with enlarged state space given by 
S = {((i,u,c), (j,v,d'))} and offspring distribution to be Poisson distributed. 
From the definition of the process d n and the relation (6.5) it is clear that 

|A«| < Y,d { r?> 9) ((u, 9 )Av,9)) < (6-8) 

U,V V 

We suppose from now on that ir has discrete support. 

Proposition 6.9 In the process d n the mean number of offspring of type (u, c") of an individual 
of type (v, c) with the parent of type (w, d) is given by 

0, if d = b, c = b, c" = g or d = g,c = b,d' = g 

m({{u'.d).{u.c)):{{u.d). (>.c"))) - <| ord = b,c = g,c" = b or d = b,c = g,c" = g 

\tt(v)(u + v), if d = b, c = b, c" = b or d = g,c = b,c" = b 

(6.10) 
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and 

m(((u',g),(u,g));((u,g),(v,b)) > ^\ir(v)v, (6.11) 
m(((u',g),(u,g));((u,g),(v,g)) < \ir(v)(u + ^) (6.12) 

Proof: For (6.10) it is enough to observe that from the definition of Color follows that is 
impossible to have green children from black parents. Since a black individual has all its children 
black, independently of its own parent type, the number of its children has the same low as 
the first generation in branching process b n defined by (6.4). The total number of children of a 
green individual has also the same low as the above one, but in this case it is possible to have 
children of both colors. Recall from (3.10) that K(b u (v)) = \-ir(v)(u + v) so that 

m(({u',g), {u,g)); ({u,g), (v,b))) +m(({u',g), (u,g)); ((u,g), (v,g))) = X7r(v)(u + v). (6.13) 

Therefore, to prove (6.11) and (6.12) it is sufficient to find the mean number of black children 
from green parent (and necessarily green grandparent): m(({u',g), (u,g)); ({u,g), (v,b))). 

Consider two incompatible rectangles R = (£,r, s,u, g) £ A2 and R' = (£' ,t' , s 1 ,u' , g) 6 Ai 
such that R G Af' , that is R is an ancestor of R'. Since R' belongs to the first generation, we 
have t' < and r' + s' > 0. By construction of the clan of ancestors A n and the branching 
process B n , rectangles R' in the first generation have Birth (R r ) < and Death (R') > 0. We 
want to compute the X(v, R, R')- number of black ancestors of f-type, namely the number 
of those V = (£v , Ty , sy , v) £ A3 with length v such that V £ Af and R £ Af'. Let 
Xb(v, R, R') and X a (v, R, R') be the number of such ancestors that died before and after time 
t = 0, respectively. 

As before, let D = D(v,R,R') = L^,u,^',u',v) x [0,r'] and 

L = L(£, u, u', v) = [max(( - u v ,£ - v), min(e + u, £' + u% (6.14) 

Then, 

X(v, R, R!) = X b (v, R, R!) + X a (v, R, R!) (6.15) 

with 

¥(X b (v,R,R') = k\u,r,T r ) = e -p»MD\ (-PbM D \) k (6 16) 
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and 



P(X a (v, R, R') = k\u, t, t') = e~^ L \ ( Pa ^ Ll) (6.17) 



where p is the probability that a rectangle who died in the area D is really an ancestral of R 
of type v, and analogously for p a . Therefore, 

Pb = n{v)F(Y - S < t) (6.18) 

= ir(v)e T {e- T ' - l)/(-r') (6.19) 

and 

p a = tt(v)¥(-S<t) (6.20) 

= 7r(w)e T (6.21) 

where Y ~ U(0, r') (death time) and S ~ exp(l) (lifetime) are independent random variables. 
We remind here that the rectangles are constructed in the negative time ("past") so r,r' < 0. 

Given R and R', X a (v, R, R') and Xb(v,R, R') are independent random variables, thus 
X(v,R,R') is Poisson distributed with mean 

n(v)(p b \\D\ + p a \\L\) = 7r(v)X\L\e^- T '\ (6.22) 

Notice that r < r' and we are given the parent relation, so r' — r has exponential distribution 
with mean one. Hence, E(e - ( r ' - ' r )) = 1/2. 
Furthermore, since £ G [£' — u, £' + u'] 

|L| =min(£ + u, £' + «')- max(£ - u) > u (6.23) 

From (6.22) and the above inequality follows 

m(((u',g),(u,g)y,((u,g),(v,b))) =E(X(v,R,R')) > \\t:(v)v (6.24) 

And we have proved (6.11) as desired. □ 

Sub-criticality: On account of the relation (6.8) we are interested in sub-criticality condi- 
tions for the green-type population of the branching process d n described above. Remember 
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that green individuals may appear in the n-th generation only as descendants of a branch made 
of green individuals only. Therefore, our aim is to establish conditions for the convergence of 
the series 

£ £ m {n \{( Wl g), (v, g)); ((v n _i , g), (u, g))) (6.25) 

n>l v n -l,u 

where m^\- ; •) = m(- ; •) is given by (6.11) and for n > 1 

mW(((«; l( /) I ( 1 ; I c));((« ll ci) I ( U /))) (6.26) 
= ^ m (n_1) (((«;,«/), (w,c)); ((u 2 ,c 2 ), (t>i, ci)))m(((v 2 , c 2 ), («i,ci)); ((ui,ci), (u,c")))- 

V2,C2 

To simplify notation let 

m*(v,«) =A7r(u)(« + S) (6.27) 
and it follows from Proposition (6.9) that for all u, t> and independently of w 

m((( w ,9),(v,g));((v,g),(u,g))) <m*(v,u). (6.28) 
Therefore, using the notation above and (3.10) we can easily estimate the series (6.25) to 

m {n) {{(w,g), (v,g)); ((v n -i,g), (u,g))) < ^ m(y, vi)m*(vi, v 2 ) ■ ■ • m*(u n _i, u) 

V„-1,U U,V1...V„-1 

= A n £ u E«! • ■ ■ 7r(«i)(« + «i)7r(«2) ( Wl + f ) ■ ■ ■ tt(u) (« n _! + |) (6.29) 

Observe that 

^ 7r(u) (v„_i + I) = + y = /* + 3i«n-i (6- 3 °) 



where /f = 4f and = 1, and define inductively 



Then we have 



7T(i; n _i + l) + — * +1 ) (/*_! + fl^i-iUn-i+l) = /* + flfan- 



(6.31) 







Pi/2 P2/2 




h 




Pi/2 p 2 /2 




Pi/2 


. 9*+i 




1 ^ . 




* 




1 ^ . 




1 



(6.32) 



and consequently, the series (6.29) is equal to 

A"2>(«l)(« + ^l)(^-l^l + = \ n (v( 9 * n - lP l + fU)+P29n-l + Pl/n-l)- 



(6.33) 
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Pi/2 P2/2 
1 Pi 



In order to find fn,9n we exponentiate T* = 
the eigenvalues of T* , e\ and £2 given by 

3/9i ± v 7 / ! + 8 /°2 



£1,2 



and two corresponding normalized eigenvectors 



£\ ~ Pi 
1 



For this operation suffices 



(6.34) 



Vi + ^i-pi) 2 

From (6.32) it follows 



V 1 + (^2-Pi) 2 



£2 - Pi 

1 



(6.35) 



1 



£1 - £2 

n* - (p n — p n 

9n — ~ ~\ £ 1 £ 2 
£l ~ £2 



£i(£i -Pi) +£%(P1 -£2)) 



(6.36) 
(6.37) 



Using the fact that |£2/£i| < 1 in (6.33) and the Cauchy-Hadamard formula, we obtain that 
the radius of convergence of the series (6.29) is 



A** 



(6.38) 



£1 3/oi + VPi + 8 P2 

Therefore, as long as A < A** the series (6.25) is absolutely convergent and consequently the 
green-type population of the process d n is almost surely finite. 

Calculation being done for countable number of types is not a real limitation. Namely, let 
V be the set of all possible types, by the same argument as used in Proposition 6.9 the mean 
number of green-type offspring in all generations is less then 



where 



m 



/ m^ l \v,du) 
n>l Jv 

f l \v,du) = \ mi" - l \v,dv')m*{v' ,du) 
Jv 



(6.39) 



(6.40) 



can be obtained inductively. 
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Suppose that the distribution of the length of the calls is absolutely continuous with respect 
to the Lebesgue measure and call ir its density. We can write 



du 

m*(v,du) = \ir(du)(— + v) (6-41) 
and the computation is completely analogous to the discrete case. 

We summarize the results proved in this section in the following theorem: 

Theorem 6.42 If A < A** , where A** is given by (6.38), then with probability one there is no 
backward oriented percolation in R. 



Remark: If n is the U(0, 1) density the condition becomes A < 1.247 
Observe the following relations: 

^->K*> , ( X ^ = IK- (6-43) 
3pi 3(pi + y/pj) 3 

The upper estimate in (6.43)is achieved in the case of just one type (fixed call length). The last 
inequality proves that a better bound for the critical value is obtained. 



7 Conclusion 

This is one of the first works where the clan of ancestors algorithm was implemented. Berthelsen 
and M0ller (2002) compared it to the dominated CFTP introduced by Kendall and Miller 
(2000). Based on simulation results, they show that the dominated CFTP is better than the 
algorithm based on the clan of ancestors in the particular case of a Strauss process 

defined on a unit square with = 100 and = (the so-called hard-core process), 0.5 
and 1 (a Poisson processes with rate 100). This is obviously the case from the description 
of the processes since the backward construction of BFA stops when the dominated Poisson 
process regenerates and usually the coupling of CFTP is achieved before it in the finite case. 
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However, it should be noticed that the algorithm based on the clan of ancestors was designed 
for sampling the infinite-volume process viewed in a finite window. This seems to be a much 
more interesting and challenging problem which has been studied in this work for the specific 
case of one-dimensional loss networks with bounded calls. No comparison was made to other 
perfect simulation schemes. 

Moreover, we can see that the simulation of the invariant measure can bring information 
about unknown variables related to the clan of ancestors. The bound described in Section 6 
was found by a better understanding of the simulation procedure. 
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